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Abstract. M87 is a nearby radio galaxy and because 
of its misaligned jet, it is possible to correlate 
detailed spatially-resolved emission regions in the 
radio, optical to X-ray waveband with unresolved 
but contemporaneous flux measurements in the TeV 
regime. Hence, M87 provides a unique opportunity 
to reveal the emission mechanisms responsible for 
high energy gamma-ray emission from active galactic 
nuclei. Observations with VERITAS since 2007 have 
resulted in 90 hours of data while 2008 observa- 
tions were part of a concerted effort involving the 
three major atmospheric Cherenkov observatories: 
H.E.S.S., MAGIC and VERITAS. As a result of the 
TeV campaign, a high flux state of M87 was detected 
in February 2008 showing multiple flares with rapid 
variability. We will present the comprehensive results 
from VERITAS observations since 2007 and also 
show preliminary results from the 2009 campaign. 

Keywords: gamma rays: observations - galaxies: 
individual (M87, VER J1230+123) 

I. Introduction 

M87 is a giant elliptical galaxy located 16Mpc away 
(redshift z = 0.00436) near the center of the Virgo 
cluster. It has been observed at all wavelengths ranging 
from radio to TeV gamma rays. Its core is an active 
galactic nucleus (AGN) powered by a supermassive 
black hole of - 3.2 x 10 9 M Q J20). The jet of M87 
does not point along our line of sight; however apparent 
superluminal motion has been observed in both radio 
lfT2l and optical [ 1 1 ] for different features along the jet, 
constraining the jet orientation to < 30° at the location 
of the knot HST-1. M87 is described as a misaligned 
BL Lac ll27ll . The proximity of M87 and its misaligned 
jet have enabled the study of its jet morphologies, 
which are similar in radio, optical, and X-rays ll22l . 
Flaring activities have been observed at these energies 
simultaneously and in different jet features [12|, which 
revealed many characteristics of relativistic jets in AGN. 

TeV emission from M87 was discovered by the 
HEGRA collaboration from their 1998-1999 observa- 
tions and was confirmed by the H.E.S.S. collabo- 
ration 1 6 1, which additionally reported year-scale and 
days-scale flux variability during a high state of gamma- 
ray activity in 2005. The observed variability timescales 
disfavor large scale gamma ray production models such 



as the dark matter annihilation model [8] and the in- 
teracting cosmic ray proton scenario [23], and favor 
the immediate vicinity of the M87 black hole as the 
TeV production site. However, the angular resolution of 
imaging atmospheric Cherenkov telescopes (IACTs) is 
insufficient to resolve any structure in M87. During the 
same period of the flare observed by H.E.S.S. in 2005, 
the Chandra X-ray observatory detected the knot HST-1 
(~ 0.8" away from the nucleus) at an intensity more 
than 50 times that observed in 2000. HST-1 was then 
suggested as a more likely source of TeV emission than 
the core lfl6l . 

The knot HST-1 has been demonstrated as a possible 
location for jet reconfinement where photons can be 
upscattered to TeV energies via the inverse-Compton 
process [25|. Several models with emission originating 
at the inner jet region have also been proposed, where 
TeV emission can be produced via inverse-Compton 
scattering [14] or via synchrotron self-Compton pro- 
cesses involving more complex jet structures |[T8l (26) • 
Leptonic models involving the electromagnetic field of 
the black hole iPTI 11241 with TeV emission coming from 
the vicinity of the black hole and not the inner jet have 
also been suggested. 

VERITAS confirmed TeV emission above 250 GeV 
from M87 in the 2007 dataset, but at a lower flux than 
what was reported by H.E.S.S. in 2005, and no vari- 
ability was detected for 2007 [3|. Day-scale variability 
was reported in 2008 by the MAGIC collaboration [7] 
with a 13-day flare during the 2008 TeV joint moni- 
toring campaign by H.E.S.S., MAGIC, and VERITAS 
0. VERITAS observation was intensified following a 
trigger alert issued by the MAGIC collaboration and 
another flare was detected [2 |. Chandra X-ray data taken 
in the same month (6-week sampling frequency) showed 
historical maximum activity coming from the core while 
the nearby knot HST-1 remained quiescent lfl5ll . The 
contemporaneous gamma-ray and X-ray flares suggest 
the core as the more probable TeV emission region, in 
contrast to the 2005 flaring activity in TeV observed by 
H.E.S.S. and in the knot HST-1 in X-rays by Chandra. 
Publication for the results from this joint campaign is 
forthcoming Q]. 

In 2009 the TeV monitoring campaign is contin- 
ued with MAGIC and VERITAS. In this paper, we 
will present preliminary results of the VERITAS 2009 
dataset, along with the results from observations begin- 
ning in 2007. 
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II. VERITAS OBSERVATIONS 

VERITAS, the Very Energetic Radiation Imaging 
Telescope Array System, is an array of four 12 m diam- 
eter imaging atmospheric Cherenkov telescopes located 
at the Fred Lawrence Whipple Observatory at Mount 
Hopkins in southern Arizona. Each telescope is equipped 
with a camera comprising 499 photomultiplier tubes 
arranged in a hexagonal lattice covering a field of view 
of 3.5°. The array is sensitive from 100 GeV to more 
than 30TeV. It has an effective area of ~ 10 5 m 2 and 
an angular resolution of ~ 0.1° (68% containment). For 
more details of VERITAS, see 0. 

M87 was observed with VERITAS for over 115 hours 
between February 2007 and April 2009 at a range of 
zenith angles from 19° to 41°. Observations in spring 
2007 were carried out during the construction phase 
and only 3-telescope data (94% of spring 2007 data) 
are used in the spectral analysis. Later observations 
(fall 2007 onward) were achieved with 4 telescopes. 
All observations were performed in wobble mode where 
M87 is tracked with a 0.5° offset to the camera center. 
After eliminating bad weather observations and unstable 
trigger rate data, over 90 hours of quality live data 
were then processed with several independent analysis 
packages |[T3l with slightly different algorithms. All 
analysis packages yield consistent results. 

Shower images are first corrected in gain and tim- 
ing using parameters obtained from the nightly laser 
calibration data. Then the images are passed through 
a two-threshold cleaning. Each shower image is then 
parametrized ifTTl . and the shower direction is recon- 
structed using the stereoscopic technique. Events are 
selected as gamma-ray like if at least two images passed 
cuts optimized for a 10% Crab Nebula flux source. The 
source region is defined by a 0.15° radius disk centered 
on the source coordinates, and all the gamma-ray like 
events within this region are summed to the ON count; 
the background is estimated from seven identically sized 
regions reflected from the source region around the 
camera center, and is summed to the OFF count iflOl . 
The ON and OFF counts are then used in the Li & 
Ma formula 17 [19| to calculate the significance of the 
excess. 

III. Results 

In 2007, M87 was detected at a statistical significance 
of 5.9 a after 44 hours of observations between February 
and April with a 3-telescope array. An average flux 
of (3.47 ± 1.12) x 10~ 12 cm _2 s _1 for energies above 
250 GeV was measured from this dataset, corresponding 
to ~ 2% of the Crab Nebula flux. No significant short- 
term flux variability was detected 0. 

In 2008, M87 was detected at 7.2 a after 41 hours 
of observations between December 2007 and May 2008. 
An average flux of (2.74±0.93) x 10 -12 cm -2 s -1 above 
250 GeV was recorded, corresponding to ~ 2% of the 
Crab Nebula flux. During 4 days in February 2008 (MJD 
54505 - 54509) flaring activity was observed in 6 hours 



of data, resulting in a 7.4 a detection and an average flux 
that corresponded to 5.3% of the Crab Nebula flux. 15 
hours of observations were performed before the flare 
period and resulted in a marginal detection (< 5er) of 
M87 at 2.0% of the Crab Nebula flux; after the flare 
period, 19 hours of observations yielded no detection of 
M87 and an upper limit of < 1.4% of the Crab Nebula 
flux. The spectral index (T) of the differential spectrum 
power-law fit with the form d^/dE = § (E/TeV)~ r 
showed no significant variation between pre-flare and 
flare period (see Table H). 

In 2009, 18 hours of observations between January 
and April yielded a marginal detection of M87 and an 
upper limit of < 1.9% Crab Nebula flux. No flaring 
activity was observed (see figure |2). 

TABLE I 

Differential spectrum power-law fit of the form 
d<$>/dE = 3> (E/TeV)- r 



data 


O 


r 




(10~ 13 cm- 2 s- 1 TeV- 1 ) 




2007 


7.4 ± 1.3 


2.31 ± 0.17 


2008-all 


5.2 ±0.9 


2.49 ±0.19 


2008-flare 


15.9 ±2.9 


2.40 ±0.21 


2008-preflare 


5.6 ± 1.5 


2.49 ± 0.26 



VERITAS 2008 
VERITAS 2007 
MAGIC 2008 
HESS 2005 
HESS 2004 




Energy (TeV) 

Fig. 1. VERITAS M87 energy spectrum of the entire dataset, in 
comparison to TeV energy spectra reported in the past. The spectral 
index of all datasets ranged from 2.22 to 2.62, and are compatible 
within statistical errors. 



IV. Discussion 

VERITAS observations of M87 spanning three ob- 
serving seasons (from 2007 to present) have shown M87 
in steady emission state and flaring state. The spectra 
obtained from these observations show no significant 
changes in the spectral index. During the 2008 joint 
monitoring campaign of M87, VERITAS observed a 
gamma-ray flare in February 2008 which spanned 4 
days, constraining the emission region size to R < 



FL, 



= ScAt/(l + z) = <5 10 



16, 



ll.UR, 



where S is the relativistic Doppler factor and R s the 
Schwarzschild radius of the M87 black hole. Rapid 
variability reported previously Q constrained the 
size of the TeV emission region to < 2.6SR S . 
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Fig. 2. Upper panel: Chandra X-ray lightcurves measured from the core and the knot HST-1. Lower panel: Nightly fluxes for energies above 
250 GeV from 2007 to present. Grey area highlights the 2008 flare period observed by VERITAS and the corresponding X-ray fluxes observed 
by Chandra. For details on the 2008 flare period, see [T] 



Even though the gamma-ray observing technique can- 
not resolve individual features of M87, the TeV emission 
size constraint has narrowed down the most probable 
TeV emission location to the unresolved core region and 
the knot HST-1. The 2008 gamma-ray flare coincided 
with the Chandra observation of historically high flux 
coming from the core while the nearby knot HST-1 
appeared to be inactive (figure |2} [1J- The contempo- 
raneous gamma-ray and X-ray flares suggest the core is 
more likely the TeV emission region, in contrast to the 
2005 flaring activity in TeV observed by H.E.S.S. and 
in the knot HST-1 in X-rays by Chandra 1T61 . 

The 2008 multi-wavelength observations of M87 in- 
cluded concurrent radio, X-ray, and TeV gamma-ray 
coverage of M87 flaring activity from the core region 
in early 2008 ||28| CD. From the 2008 multiwavelength 
data, the TeV emission region is likely the unresolved 
core. However, the knot HST-1 is still a possible can- 
didate for the 2005 flare. Current models do not favor 
one over the other, and both the core and HST-1 remain 
as candidates for TeV emission. As of the end of April 
2009, the 2009 monitoring work has shown no flaring 



activity from M87. Further multi-wavelength monitoring 
can potentially provide additional constraints on the 
environment of M87 and more insights into the emission 
mechanism of AGN. 
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